Toroidal normal forms for bifurcations in retarded 
functional differential equations I: Multiple Hopf and 
transcritical/multiple Hopf interaction 

Younsim Choi 
Department of Mathematics and Statistics 
University of Ottawa 
Ottawa, ON KIN 6N5 
CANADA 

Victor G. LeBlanc 
Department of Mathematics and Statistics 
University of Ottawa 
Ottawa, ON KIN 6N5 
CANADA 

May 18, 2005 



Abstract 

For finite-dimensional bifurcation problems, it is well-known that it is possible to 
compute normal forms which possess nice symmetry properties. Oftentimes, these sym- 
metries may allow for a partial decoupling of the normal form into a so-called "radial" 
part and an "angular" part. Analysis of the radial part usually gives an enormous 
amount of valuable information about the bifurcation and its unfoldings. In this pa- 
per, we are interested in the case where such bifurcations occur in retarded functional 
differential equations, and we revisit the realizability and restrictions problem for the 
class of radial equations by nonlinear delay-differential equations. Our analysis allows 
us to recover and considerably generalize recent results by Faria and Magalhaes [11] 
and by Buono and Belair [4]. 
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1 Introduction 



Delay-differential equations are used extensively in the modeling of a multitude of phenomena 
in the life sciences [3, 20, 22], physics [19, 27], atmospheric sciences [25], engineering [24], 
economics [28] and beyond. This has motivated a flurry of activity on the mathematical 
side to try to understand the behavior of this class of equations and to develop a theoretical 
framework suitable for their analysis. 

It is now well-understood that retarded functional differential equations (RFDEs), a class 
which contains delay-differential equations, behave for the most part like infinite-dimensional 
ordinary differential equations. The upshot is that many of the techniques and theoretical 
results of finite-dimensional geometrical dynamical systems are portable to RFDEs. In par- 
ticular, versions of the stable/unstable and center manifold theorems in neighborhoods of an 
equilibrium point exist for RFDEs [18]. For example, near a bifurcation point in a RFDE, 
the flow is essentially governed by a vector field on an invariant center manifold. This has 
allowed for the successful application of the vast machinery of bifurcation theory to many 
problems which are modeled by RFDEs, e.g. [2, 24]. Parallel to this, techniques for sim- 
plifying vector fields via normal form changes of coordinates have been adapted to RFDEs 
[8, 9], and has allowed for further insight into the qualitative behavior of RFDEs. 

This paper is concerned with the bifurcation theory of RFDEs. In particular, we will be 
interested in the so-called realizability problem for normal forms of vector fields which arise 
via center manifold reduction of RFDEs. 

Realizability problem: 

Suppose B is an arbitrary mxm matrix. For the sake of simplicity, suppose additionally that 
all eigenvalues of B are simple. Let C([—r, 0], R) designate the space of continuous functions 
from the interval [— r, 0] into R, and for any continuous function z, define z t G C([—r, 0], R) 
as z t (9) = z(t+9), —r<9<0. It is then possible [10] to construct a bounded linear operator 
L : C([— r, 0], R) — > R such that the infinitesimal generator A for the flow associated with 
the functional differential equation 

z{t)=Lz t (1.1) 

has a spectrum which contains the eigenvalues of B as a subset. Thus, there exists an Tri- 
dimensional subspace P of C([—r, 0], R) which is invariant for the flow generated by A , and 
the flow on P is given by the linear ordinary differential equation (ODE) 

x = Bx. 

In our case, we will be especially interested in the case where the eigenvalues of B all have 
zero real parts, and the spectrum of A does not contain any elements with zero real part 
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other than those which belong to the spectrum of B. 

Now, suppose (1.1) is modified by the addition of a nonlinear delayed term 

z(t) = L z t + az(t + t) 2 , (1.2) 

where a G K is some coefficient and r G [— r, 0] is the delay time. Then the center manifold 
theorem for RFDEs[18] can be used to show that the flow for (1.2) admits an m-dimensional 
locally invariant center manifold on which the dynamics associated with (1.2) are given by 
a vector field which, to quadratic order, is of the form 

x = Bx + ag(x), (1.3) 

where g : W 71 — > R m is a fixed homogeneous quadratic polynomial which is completely 
determined by L and r, and a is the same coefficient which appears in (1.2). We immediately 
notice that for fixed L and r, (1.3) has at most one degree of freedom in the quadratic term, 
corresponding to the one degree of freedom in the quadratic term in (1.2). However, whereas 
one degree of freedom is sufficient to describe the general scalar quadratic term involving 
one delay in (1.2), it is largely insufficient (if m > 1) to describe the general homogeneous 
quadratic polynomial / : lR m — > R m . Therefore, there exist m-dimensional vector fields 
x = Bx + f(x) (where / is homogeneous quadratic) which can not be realized by center 
manifold reduction (1.3) of any RFDE of the form (1.2). One quickly notices that the 
situation could be improved if we allow the nonlinear terms in (1.2) to depend on more than 
one delayed times, i.e. 

2 

z(t) = L(z t ) + «Ui*~iM l + r i)) n • • • (*(* + T j)Y j > (1-4) 

iiH Mj=2 

where the a^...^. are real coefficients and Ti, . . . , tj G [— r, 0] are the delay times. The center 
manifold equations for (1.4) truncated to quadratic order are 

2 

x = Bx+ ^ a ili2 ...i j g ili2 ...i j {x), (1.5) 

il,...,ij= 

iiH hij=2 

where gi^...^ : M m — > M. m are fixed homogeneous quadratic polynomials which are com- 
pletely determined by L and Ti,...,Tj. Thus, the subspace of m-dimensional vector fields 
x = Bx + f(x) (where / is a homogeneous quadratic) of the form (1.5) is potentially larger 
than those of the form (1.3). Of course, there is nothing particularly special about the 
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quadratic order, and one could repeat the above discussion to include progressively higher 
order nonlinearities. Without loss of generality, we could also limit our attention to only 
those / which are in normal form with respect to the matrix B. The particular version of 
the realizability problem which will interest us in this paper is the following: 

Given: 

• anmxm matrix B whose spectrum consists solely of simple eigenvalues with zero real 
parts, 

• a bounded linear operator L : C([— r, 0],R) — > R such that the infinitesimal generator 
A for the flow associated with the functional differential equation (1.1) has a spectrum 
which contains the eigenvalues B as a subset, and no other part of its spectrum on the 
imaginary axis 

• an integer t > 2 

• a polynomial f : M. m — ► M. m of degree £ such that /(0) = and Df(0) = 0, and f is 
in normal form with respect to the matrix B 

does there exist an RFDE of the form 

z(t) = Lz t + F(z(t + T 1 ),...,z{t + T j )), (1.6) 

such that the center manifold equations for (1.6), in normal form and truncated to order t, 
are x = Bx + f{x) ? 

This question was answered in the affirmative for scalar RFDEs [16, 17] and in [10] for 
n-dimensional RFDEs in general. In the scalar case (which will be of interest to us), the 
result states that there is generically a solution to the realizability problem for general / as 
stated above if j (the number of distinct delays in (1.6)) is at least equal to m (the dimension 
of the center subspace). 

The main purpose of this paper is related to the optimality of the above sufficient number 
(j = m) of delays, in light of some recent results by Faria and Magalhaes [11], and by Buono 
and Belair [4], which we now describe. 

Simple Hop f, (0,±iu), and (± iu>i, ± ilo-i) bifurcations: 

In [11] and [4], the authors consider the optimality of the solution j — m to the realizability 
problem for scalar RFDEs in some special cases. Consider one of the following three separate 
cases for the matrix B: 
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• B is a 2 x 2 matrix whose eigenvalues are ±iu, u > 0, 

• 5 is a 3 x 3 matrix whose eigenvalues are and ±iu, u > 0, 

• B is a 4 x 4 matrix whose eigenvalues are ±iu>i, ±icj 2 , where lo\ > and cj 2 > are 
rationally incommensurate. 

Let L : C([—r, 0], R) — > R be a bounded linear operator such that the infinitesimal generator 
A for the flow associated with the functional differential equation (1.1) has a spectrum 
which contains the eigenvalues of B as a subset, and has no other part of its spectrum on 
the imaginary axis. Therefore, a general RFDE of the form 

z(t) = Lz t + N(z t ), (1.7) 

where N(0) = 0, DN(0) = has an equilibrium solution z = undergoing respectively 
a simple Hopf bifurcation, a steady-state/Hopf interaction, or a non-resonant double Hopf 
bifurcation. Normal forms and versal unfoldings for each of these bifurcations are well- 
known. In the first case, using normal form changes of coordinates and then converting to 
polar coordinates X\ = pcos 0, x 2 = psin for the center manifold, the normal form (to 
cubic order) is 

p = ap 3 , = u + bp 2 . (1.8) 

If the coefficient a in (1.8) is non-zero, then the higher-order terms have no qualitative 
effects. In this case, the p equation completely determines the bifurcation. Now, from the 
above-mentioned solution j — 2 to the realizability problem [16, 17], we can conclude that 
if N(zt) is of the form N(z t ) = F(z(t + Ti), z(t + r 2 )), then any value of a and b in (1.8) can 
be realized by means of center manifold reduction of (1.7). However, Faria and Magalhaes 
show that, in fact, any value of a in the determining p equation of (1.8) can be generically 
realized if N(z t ) only involves one delay, i.e. N{z t ) = F(z(t + Ti)). Similarly, they show that 
the versal unfolding 

p = \p + ap 3 

of the p equation in (1.8) can be realized by a RFDE of the form z(t) = L(z t ) + A z(t + T\) + 

F(z(t + n)). 

For the (0, ± iuo) case, the center manifold is three-dimensional. Normal form changes of 
coordinates and the use of center manifold coordinates x±, x 2 = pcos 0, x 3 = psin yield 
the following equations (to quadratic order) 

xi = bix\ + b 2 p 2 

p = a 1 x 1 p 

= u + 0(\ Xl ,p\ 2 ). 
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If the coefficients bi, b 2 and a\ satisfy certain generic non- degeneracy conditions, then the 
higher-order terms have no qualitative effects. The bifurcation is thus characterized by the 
ii, p subsystem. In this case, Faria and Magalhaes show that any value of b\, b 2 and a\ can 
be realized by center manifold reduction of an RFDE involving only 2 delays (which is less 
than the predicted value (j — m — 3) from the solution to the realizability problem [16, 17]). 

Finally, a similar result was shown for the non- resonant double Hopf bifurcation in [4]: 
whereas the solution of the realizability problem [16, 17] predicts that j = m = 4 delays are 
sufficient, it is shown in [4] that 2 delays are sufficient to realize, to cubic order, the "radial 
part" (p'i, p 2 ) of the center manifold equations 

Pi = (a*i + anPi + a 12 (%)pi 

Pi = (P2 + a 2 ip\ + a 22 p\)p 2 

6>i = io\ + 0(|/ii,// 2 |, |pi,P2| 2 ) 

O2 = U 2 + 0(\p 1 ,p 2 \, \P1,P2\ 2 )- 

In all three cases above, it is also shown that this smaller number of delays (1 in the 
simple Hopf case and 2 in both the (0, ±icu) and (± iu>i, ± icu 2 ) cases) is optimal, in the 
sense that anything less will lead to restrictions on realizability of the various coefficients 
which appear in these normal forms, and consequently to restrictions on the possible phase 
portraits in the classification of the versal unfolding of these respective singularities. 

Overview: 

The questions of realizability and restriction for normal forms and unfoldings of bifurcations 
in RFDEs are particularly important from a modeling point of view. Indeed, given a specific 
RFDE model (perhaps depending on many parameters) undergoing a local bifurcation, it is 
important to be able to characterize the range of possible dynamics accessible from within the 
model near the bifurcation point. From our discussion above, we see that knowledge of the 
abstract finite-dimensional bifurcation problem and its unfoldings is not sufficient in general 
to answer this question. Indeed, the specific form of the RFDE may restrict this range of 
possible dynamics. Depending on the functional form of the RFDE (e.g. how many distinct 
delays are involved), some phase diagrams which are possible in the unfolding of this given 
bifurcation may not be realizable in the RFDE. This could have important consequences in 
the interpretation of the model, especially as it pertains to the actual phenomenon being 
modeled. 

The purpose of this paper is to further study these issues of realizability and restrictions, 
in light of the previously discussed results in [11] and [4]. We will develop a unified theoretical 
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framework for these results, which will consequently allow for considerable generalizations 
of these results. 

Specifically, we will exploit and generalize the following common elements of the three 
specific cases studied in [11] and [4]: it is possible to make a canonical choice of normal 
form transformations on the center manifold equations which lead to a normal form with 
nice symmetry properties - it is equivariant with respect to an action of a torus group. This 
toroidal equivariance can be used to achieve a partial decoupling of the normal form into a 
"radial part" and an "angular part" . In many cases (in particular, in the three specific cases 
studied in [11] and [4]), the radial part characterizes the essential features of the bifurcation. 
It is then reasonable to 

investigate the realizability problem for the radial part of the normal form and its unfoldings, 

which is the goal we seek. Along with developing a theoretical framework to achieve this 
goal, we will in fact make the following generalizations to the results of [11] and [4]: 

• we will assume that the spectrum of the matrix B consists of simple eigenvalues, and 
has one of the two following forms 

spec(.B) = {± icui, . . . , ± icUp} or spec(B) — {0, ± iui, . . . , ± icu p }, (1.9) 

where u>i, . . . , cu p > are independent over the rationals, 

• we will investigate realizability of the radial part of the normal form to any order, and 
not just quadratic or cubic. This is important in cases where nonlinear degeneracies 
are present. 

Similarly to [4, 11], we will limit our analysis to the case of scalar RFDEs. While studying 
the realizability problem in the context of general n > 1 dimensional systems of RFDEs is 
certainly very important, our computations indicate that there is an enormous increase in 
algebraic complexity involved. Consequently, a unified concise framework allowing for the 
simultaneous analysis of all cases (scalar and systems) appears at this point to be a difficult, 
albeit not impossible, goal to achieve. 

In the second case for the spectrum of the matrix B in (1.9), there is an technical subtlety 
which arises as it pertains to the possible unfoldings of this singularity. In fact, there are two 
algebraically different ways to construct an unfolding, depending on whether the "steady- 
state" mode (correspongind to the eigenvalue) in the interaction is of saddle-node type 
or of transcritical type. It turns out that the transcritical case can be treated in the same 
framework as the first case for spec(-B) in (1.9), but it would be extremely cumbersome 
to attempt to treat the saddle-node case within this same framework. Therefore, we have 
chosen to treat the saddle-node case in a separate paper [5]. 
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This paper is organized as follows. In Section 2, we will give a brief review of the theory 
of the center manifold reduction and normal form transformations of RFDEs as developed 
by Faria and Magalhaes in [8, 9]. In Section 3, we will review how to make a canonical 
choice of normal form which possesses useful toroidal symmetry properties, and exploit this 
symmetry to achieve a partial decoupling of the normal form into a "radial part" and an 
"angular part". The radial part possesses residual reflectional-type symmetry, which will 
be important for the subsequent analysis. In Section 4, we set a framework and establish 
an important surjectivity result which will be crucial to study the realizability problem for 
the radial part. Our main results on realizability and restrictions are given in Section 5. In 
this section, we also give some results which hint at the rudiments of a singularity theory 
for RFDEs. In Section 6 we show how the specific results of [11] and [4] are recovered by 
our main results. We end with some concluding remarks in Section 7. Some of the proofs in 
Section 3 are relegated to the appendices. 

2 Functional Analytic Framework 

In this section we will briefly recall some standard results and terminology in the bifurcation 
theory of RFDEs in order to establish the notation. For more details, see [8, 9, 18]. 

2.1 Infinite dimensional parameterized ODE 

Suppose r > is a given real number and C = C ([— r, 0] , R) is the Banach space of 
continuous functions from [— r, 0] into R with supremum norm. We define Zt G C as 
Zt (9) — z(t + 0) , — r < 9 < 0. Let us consider the following parameterized family of nonlin- 
ear retarded functional differential equations 

z{t)=L{n)zt + F{zt,n), (2.1) 

where L : C x R s — > R is a parameterized family of bounded linear operators from C into R 
and F is a smooth function from Cxi s into R. In this paper, we will assume the following 
hypothesis on F 

Hypothesis 2.1 F(0, 0) = and DF(0, 0) = 0. 

A consequence of Hypothesis 2.1 is that in a Taylor expansion of (2.1), there are no terms 
which are Zt independent and linear in //. While this is not a restriction in one of the cases 
we will be studying in this paper (multiple non- resonant Hopf bifurcation), it is a restriction 
in the other case (steady-state/multiple non- resonant Hopf interaction). Note however that 
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Hypothesis 2.1 includes as a special case the physically interesting case in which z = is 
an equilibrium for (2.1) for all /i, which is the case of interaction between a transcritical 
bifurcation and multiple non-resonant Hopf bifurcation (see [21]). In a sequel to this paper, 
we will relax Hypothesis 2.1 to simply F(0,0) = and DiF(0, 0) = 0, which is the generic 
saddle-node case. This relaxation of Hypothesis 2.1 leads to some technical complications 
which would make a unified treatment of both cases simultaneously extremely cumbersome 
and lengthy. Therefore, for the sake of clarity, we have decided to treat these two cases 
separately (see [5]). 

The bounded linear map L(fi) can be represented in an integral form as 

L(/i)0= f [dr)„ (9)} (0) , 

J — r 

where r/^ (9) is a measurable function on [— r, 0]. Denote L = L(0), and rewrite (2.1) as 

z (t) = L z t + (L(/i) - L )zt + F(z t , fi) = L z t + F(z t , /i), (2.2) 

where F(z t ,fi) = (L(/i) - L )z t + F(z t ,fi). 

Let A(n) be the infinitesimal generator of the flow for the linear system z = L(fi)zt, with 
spectrum <j(A(fi)), and denote by A M the set of eigenvalues of a(A(fj,)) with zero real part. 

The set A , which consists of the roots of the characteristic equation 

o 

det A (z) = 0, A (z) = z - J [dr] (9)] e zd , (2.3) 

— r 

with zero real part will play an important role. 

Hypothesis 2.2 Throughout the rest of the paper, we assume the following hypotheses on 
A M and Aq : 

(a) Card(A^) < Card(A ) for ji small, 

(b) Each element of A is a simple eigenvalue of A(0), and A has one of the following 
two forms: 

A = {±iuii, . . . , ±iu p } (multiple non-resonant Hopf bifurcation), or 

A = {0, ±iui, . . . , ±iu p } (steady-state/multiple non-resonant Hopf interaction), 

where ui,...,u p , are independent over the rationals, i.e. if r 1; . . . ,r p are rational 
numbers such that Y7j=i r j UJ j = ®> then r\ — • ■ ■ — r p — 0. 
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Let P be the invariant subspace for A = A(0) associated with the eigenvalues in A , 
and let $ = (ipi . . . ip m ) be a matrix whose columns form a basis for P. 

In a similar manner, we can define an invariant space, P*, to be the generalized eigenspace 
of the transposed system, Aq associated with A , having as basis the rows of the matrix 
\J/ =col(/0i, . . . , tpm). Note that the transposed system, Aq is defined over a dual space 
C* = C ([0, r] , R) , and each element of * is included in C*. The bilinear form between C* 
and C is defined as 

o e 

(v, <t>) = v (o) 4> (o) - y yw - *) [**> (*)] (o dc (2.4) 

-r 

Note that $ and \l/ satisfy $ = 5$, ^ = — ^B, where flisanmxm matrix whose spectrum 
coincides with A . 

We can normalize ^ such that (\l/, $) = J, and we can decompose the space C using the 
splitting C = P © Q, where the complimentary space Q is also invariant for A . 

Faria and Magalhaes [8, 9] show that (2.1) can be written as an infinite dimensional 
ordinary differential equation on the Banach space BC of functions from [— r, 0] into R which 
are uniformly continuous on [— r, 0) and with a jump discontinuity at 0, using a procedure 
that we will now outline. Define X to be the function 



X o (0) 



9 = 

-r < 6 < 0, 



then the elements of BC can be written as £ = ip + X \, with <p e C and A e R, so that 
is identified with C x R. 
Let 7r : .BC — > P denote the projection 

n(<p + X \) <p) + *(0)\], 

where (f E C and A G R. We now decompose z t in (2.1) according to the splitting 

BC = P©ker7T, 

with the property that Q C ker7r, and get the following infinite-dimensional ODE system 
which is equivalent to (2.1): 

x = Bx + V(0)[(L(n)-L o )(<Z>x + y) + F(<Z>x + y,n)] 

d ( 2 - 5 ) 
— y = A Q iy + (I-n)X [(L( t i)-L )($x + y) + F(<S>x + y, f i)], 
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where x G R m , y G Q l = Q H C 1 , (C 1 is the subset of C consisting of continuously differen- 
tiable functions), and Aqi is the operator from Q 1 into ker7r defined by 

A Q np = if + X [Lif - (p 



2.2 Faria and Magalhaes normal form 

Consider the formal Taylor expansion of the nonlinear terms F in (2.2) 



i>2 



G R s , 



where Fj(w) = Hj(w, . . . ,w), with Hj belonging to the space of continuous multilinear 
symmetric maps from (C x R s ) x • • • x (C x R s ) (j times) to R. If we denote fj = (fj, ff), 
where _ 

fj(x,y,fi) = VMFjiQx + y,^ 

ff(x,y,n) = (I-ir)X Fj(^x + y,fj,), 
then (2.5) can be written as 

x = Bx + Y fj(x,y,^) 

d (2-6) 
jV = A Q iy + 22fj(x,y,n) 

i>2 

The spectral hypotheses we have specified in Hypothesis 2.2 are sufficient to conclude 
that the non-resonance condition of Faria and Magalhaes [8, 9] holds. Consequently, using 
successively at each order j a near identity change of variables of the form 

(x,y) = (x,y) + U j (x, f i) = (x, y) + (Uj(x, //), Uf(x, //)), (2.7) 

(where Uj' 2 are homogeneous degree j polynomials in the indicated variables, with coefficients 
respectively in R m and Q l ) system (2.6) can be put into formal normal form 

x = + gj(x,y,n) 

d (2-8) 
j t V = A Q iy + 2_^g](x,y,fi) 

i>2 
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such that the center manifold is locally given by y = and the local flow of (2.1) on this 
center manifold is given by 



i>2 

The nonlinear terms in (2.9) are in normal form in the classical sense with respect to the 
matrix B. 

3 Bifurcations with Toroidal Normal Forms 

With equation (2.9) in mind, in this section we will discuss normal form transformations of 
the general parameterized system 



where the spectrum A of the matrix B is as in Hypothesis 2.2(b). As much as possible, we 
will treat both cases of Hypothesis 2.2(b) (i.e. whether or not A includes 0) simultaneously 
by adopting a notation which uses integers k and d, which should be interpreted as having 
the values k = 2p and d — p in the case where A = . . . , ±ioo p }, and the values 

k = 2p + 1 and d — p + 1 in the case where A = {0, ±iu>i, . . . , ±iu p }. 

It will be extremely useful to use complex coordinates for the last 2p components of the 
space R K , so that we can identify 




(2.9) 



x 



Bx + f(x,fi) 



(3.1) 



fi = 0, 



{ (X U X U . . . ,X p ,X p ) I Xj G C, j = 1, . . .,p} 



if k = 2p 



K 



{ (x ,x 1 ,x 1 , . . . ,Xp,x p ) I x e R, Xj e C, j = 1, . . . ,p} if n = 2p+l. 



Then, without loss of generality, we may assume that 



B = diag(ia>i, —ioui, . . . , iu p , —iu p ) or B = diag(0, iu>i, iuii, . . . , iu p , —iu p ) 



(3.2) 



depending on which case of Hypothesis 2.2(b) is being considered. 
At times, it will be convenient to write (3.1) as 



x 



Bx + f(x), 



(3.3) 





(3.4) 
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The section is divided into four subsections. In the first, we will give a brief review of 
results on symmetric normal forms with parameters. Most (if not all) of these results are 
largely well-known in the unparameterized case (see for example [7, 13]), and only minor 
modifications are required to obtain the parameterized versions we present herein. 

In the second subsection, we will define an equivariant projection operator which will be 
useful in the computation of symmetric normal forms. 

In the third subsection, we will specify how the symmetry of these normal forms can be 
exploited in order to achieve a partial decoupling of the normal form. 

Finally, in the fourth subsection, we will introduce a splitting of our spaces of polynomials 
which naturally decomposes any vector field into a singular parameter independent part plus 
a perturbation. 

3.1 Normal forms and toroidal symmetry 

For B as in (3.4), let B l denote the transpose of B and let Y = {e sBt \ s G R} (where the 
closure is taken in the space of (k + s) x (k + s) matrices), and note that T is an abelian 
connected Lie group isomorphic to TP, where T p is the p-torus: 

f {diag(e ifll ,e- ifll ,...,e i0 p,e- ifl p,l,...,l) | 6j G S 1 , j = 1, . . . ,p } if re = 2p 
T p = < 

{ {diag(l,e if?1 ,e- iei ,...,e^,e-^,l...,l) | 8 3 r e S 1 , j = 1, . . . ,p } if n = 2p+l. 

(3.5) 

Definition 3.1 For a given integer £ > 2, a given normed space X, and for k — 2p (respec- 
tively k = 2p + 1), we denote by H£ +S (X) the linear space of homogeneous polynomials of 
degree £ in the k + s variables x = (x±,Xi, . . . , x p , x£) (respectively x = (x , Xi,x~[, . . . , x p , x£) ) 
and fi = (fii, . . . ,/j, s ) with coefficients in X. For X = R K+S , define H^ +S (R K+S ,T) C 
H^ +S (M, K+S ) to be the subspace of T -equivariant maps, i.e. 

f E H* +S (R K+S ,T) ^ 

f G H* +S (R K+S ) and ifin^Z) = f(x), Vx = (x,fi) G R K+S , V 7 G T. 

For the general class of near-identity changes of variables + h{x) for (3.3), it is well- 

known that we can eliminate from (3.3) all nonlinear terms which are in the range of the 
homological operator 

Cb : H* +S {R K+S ) — ► H^ +S (R K+S ) 

„ „ (3-6) 
/ ^ (jC § f)(x) = Df{x)Bx - Bf(x). 
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Thus, we must define in H^ +S {R K+S ) a complimentary space to range £g. Of course, such a 
space is not unique. However, there exists a nice canonical choice which will be extremely 
useful for our purposes (see for example [7, 13]). 



Proposition 3.2 



H, 



k+s fjaK+s 



f ,r) © range £j 



The usefulness of Proposition 3.2 is that it is straightforward to compute the general element 
of H* +S (R K+S ,T). 

Lemma 3.3 Let B be as in (3.3). Then a smooth vector field f : R K+S — ► R K+S is T p - 
equivariant if and only if f has one of the following forms 



( a 1 (x 1 x 1 ,...,x p x p ,fi)x 1 \ 
ai(x{xl, . . .,x p Xp~, /i) x x 



Qjp{x\X\^ • • • i $p%pi /^) Xp 

cip^XiXii • • • i XpXpi /x) Xp, 
bi(x{xi, . . .,x p Xp~,ii) 



\ bs(x\Xi, . . . , XpXp, /i), J 



or 



/ ao(xo,XiXi, . . . ,x p Xp,fx) \ 
ai(x , x{xi, . . . , XpXp~, fi) xi 

di{Xo, X±Xi, . . . , XpXp, /i) X\ 



QpyKQi •''l^'l) • • • j •Kp-Kpi AO -Epi 

6i(x ,xixT, . . .,x p Xp~,fx) 



\ b s (x ,x 1 xi,...,x p x p ,fj), J 



(3.7) 



respectively if k — 2p or k = 2p + 1, where ai, . . . , a p are smooth and complex-valued, and 
ao, bi, . . . b s are smooth and real-valued. 

Proof This is a standard result which is a consequence of Schwarz lemma [23]. See also 
[13]. ■ 

Proposition 3.2 and Lemma 3.3 are not exactly in a form suitable for our purposes, since 
the vector field / in (3.3) has the special form / = (/, 0) which we require our normal form 
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changes of variables to preserve. Since we are only interested in the first re components of 
(3.3), we would like to obtain a splitting of H* +S (R K ) akin to the splitting of H£ +S (R K+S ) in 
Proposition 3.2. For this purpose, we will need the following 

Definition 3.4 



(a) We define H* +S (R K , T p ) to be the subset ofH* +s (R K ) consisting of mappings f : R K+S — 
whose components are of the form of the first re components of (3.7). Note that 
H^ +S (R K , T p ) consists precisely of the T p -equivariant elements of H£ +S (R K ); that is, 

f e H* +S (R K ,T P ) 

/ E H« +S (R K ) and /( 7 o^) = lof(x,fi), V 7 o G T , V(x,/x) G R K+S , 

where T is the group of rex re matrices which is isomorphic to T p , and is parameterized 
as 

{ diag(e if \ e" if \ . . . , e i0 ^ , e~ i0 ") \ Bj G S 1 , j = 1, . . .,p} if re = 2p 

{diag(l,e iei ,e-^,...,e^,e-^) | Oj G S 1 , j = 1, . . . ,p } if re = 2p+l, 

(3.8) 

(b) We define the following operator 

C B : H^ S {R«) — > ij;+ s (M«) 
/ (C B )(f)(x,ii)=D x f(x,ii)Bx-Bf(x,ii). [6 ^ ) 

iVote too* £s(/,0) = (£ B /,0). 
Proposition 3.5 

#; +s (M K ) = #; +s (M K , T p ) © range £ B . 



Proof The proof is given in the appendix. 



3.2 Equivariant projection 

In this section, we will construct an appropriate linear projection associated with the splitting 
of H^ +S (R K ) given in Proposition 3.5. This projection has very nice algebraic properties, and 
will be useful when we prove our main results later. 
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Definition 3.6 Let / denote the normalized Haar integral on T = T p (see (3.8)). We 
define the linear operator 

A : H^ +S (R K ) — ► H* +S (R K ) 
f^(Af)(x,»)= [ 1 f{ 1 ~ 1 x^)d 1 

J To 

Proposition 3.7 A is a projection. Furthermore, 

ranged = H« +S (R K , T p ) (3.10) 

and 

ker A = range Cb (3-11) 

Proof The proof is given in the appendix. 

Since A is a projection, then H^ +S (R K ) = ker A © ranged, and Proposition 3.7 shows 
that this decomposition is precisely the decomposition of H^ +S (R K ) given in Proposition 3.5. 
Thus, for any / G H* +S (R K ), write 

f = Af + (I-A)f, 



and note that Af is T p -equivariant and that (/ — 
exists h G H* +S (R K ) such that C B h = (I - A)f. 

3.3 Phase decoupling 



A)f G ker A. From Proposition 3.7, there 



The following example serves as an illustration of a trivial (well-known) case in which normal 
form toroidal symmetry leads to a decoupling of the equations in the normal form. 

Example 3.8 In the case where B = diag(icj, — iuo), the normal form has the rotational 
symmetry of a one- dimensional torus: ( '), 9 G T 1 , and it is easy to verify 

that the most general T 1 -equivariant differential equation has the form 

x = f(xx)x, (x G C) (3.12) 

and its complex conjugate where f is complex-valued, and /(0) = ioj. Writing x = re l9 leads 
to the equations 

f = Re(/(r 2 )) r, 6 = Im(/(r 2 )). (3.13) 
16 



We note that 6 does not appear in the r equation, and that the r equation has a reflectional 
symmetry r — > —r. The analysis of the normal form (3.12) then essentially reduces to a one- 
dimensional problem (the f equation in (3.13)) which possesses some residual (reflectional) 
symmetry. 

In fact, this example is a special case of a more general result which holds when the 
spectrum of B satisfies Hypothesis 2.2(b), and which we now outline. 
From Lemma 3.3, we get the following 

Corollary 3.9 Suppose the spectrum of B satisfies Hypothesis 2.2(b). Then a smooth vector 
field f : R K+S — ► 1R K with /(0, 0) = 0, Df(0, 0) = is T = JP-equwanant if and only if f 
has the form 

( ai (x{x~[, . . . , x p Xp~, /i) xi \ 



ai{xix 1} . . . , x p x p , fx) xi 



if k = 2p 



a,p{x\X\^ • • • , XpXp, /i) Xp 

y Clp^XiXi, . . . , XpXp, /i) Xp, J 

f(x,fi) = { ( a (x ,x 1 xJ,...,x p Xp',ii) \ (3.14) 
ai{Xo, x±x±, . . . , XpXp, n) X\ 

0\{x,q, X\X\, . . . , XpXp, /i) X\ 

if K = 2p+1, 

Q'pi^O'} X\X\, • • • , XpXp, ^x) Xp 
k \ Qip\%Q, X\X\, • • • , XpXp, /i) Xp, J 

where a\, . . . , a p are smooth and complex-valued, and a is smooth and real-valued. 

Proposition 3.10 Consider a differential equation x = Bx + f(x,fi), where f is as in 
(3.14). Then under the under the change of variables Xq = po, Xj = Pje l9j , j = 1, . . . ,p, this 
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differential equation transforms into 

Pj = ~Re(aj{p\, . . . , p 2 p , n)) pj, j = l,...,p if « = 2p 



Po = ao(po, (%,...,($, fi) (3.15) 

if K = 2p+1, 

Pj = Re(aj(p ,pl,...,pl,p)) pj, j = l,...,p 
and 

{lm(aj(pl, . . . , p 2 , p)), j = l,...,p if n = 2p 
(3.16) 
lm(aj(po,pl,...,pl,p)), j = l,...,p if K = 2p+1. 

Proof This is a simple computation. ■ 

We will call the subsystem (3.15) the uncoupled radial part of the normal form (3.14). For 
many practical purposes of interest, it is sufficient to consider only the uncoupled radial part 
(3.15) in the analysis of (3.14). For example, small-amplitude equilibria of (3.15) correspond 
to periodic solutions or invariant tori of the full normal form (3.14). Oftentimes, given some 
normal hyperbolicity conditions, these invariant objects for (3.14) persist as invariant objects 
in the original system (3.1). In fact, in the case of non- resonant double Hopf bifurcation (A = 
{±iuJi, ±iuj 2 }) and in the case of saddle-node/Hopf interaction (A = {0, ±iu 1 }), it is well- 
known [26] that given some generic non- degeneracy conditions on the coefficients of the lower- 
order nonlinearities, the radial equations (3.15) (suitably truncated) completely determine 
the dynamics in the full system (3.1) up to topological equivalence. So, it is reasonable to 
investigate the realizability of the uncoupled radial part (3.15) by center manifold reduction 
(2.9) of the RFDE (2.1). 

We now introduce an integer d which should be interpreted such that d = p in the case 
where A = {±iu\, . . . , ±iu p }, and d = p + 1 in the case where A = {0, ±iui, . . . , ±iu> p }. 
Denote by Z 2jP the group whose action on M. d is given by 

(p 1 ,...,p p ) -y (Aipi,...,A p pp) if d = p 

(3.17) 

(p , P!, . . . , p p ) -> (po,Aipi,...,A p pp) if d = p + l, 
where Xj E {1, -1}, j = 1, . . . ,p. 

Definition 3.11 For a given integer £ > 2, a given normed space X, and for d = p (re- 
spectively d = p + 1), we denote by Hf +s (X) the linear space of homogeneous polynomials 
of degree £ in the d + s variables p = (pi,...,p p ) (respectively p = (po, p\, . . . , p p ) and 
p = (px, ...,p s ) with coefficients in X. Denote by Hf +s (R d ,Z 2:P ) C Hf +s (R d ) the subspace 
of Hf +s (M. d ) consisting of Z2 tP -equivariant polynomials. 
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It is easy to show (see [13]) that the most general element of H d+s (R d , Z 2)P ) has the form 

/ h 1 (p 2 ,...,p 2 p ,p)p 1 \ 
\ Kipi, ■ ■ ■ » p%p) P P , J 



if d — p 



( h (p , p 2 , . . . , p 2 p , p) \ 
h(p , pi, ■ ■ ■ , p 2 p , p) pi 

\ h P {po, pi, ■ ■ ■ , P 2 P , P) P P , J 



if d = p + 1 



and one immediately notices the similarity with (3.15). It then becomes useful to define the 
following surjective linear mapping 



n : H* +S (R K ,T P ) — ► H d+s (R d ,Z 2tP ) 



(3.18) 



which is defined by sending the general element (3.14) of H^ +S (R K ,T P ) to the following 
element of H d+S (R d , Z 2 , p ) : 

/ Re(a 1 (p\, . . . , p 2 p , p)) p x \ 

\ RcMpI, ■ ■ ■ , P 2 P , P)) P P J 

( a (p , pj, . . . , p 2 p , p) \ 
Re(ai(p ,P?,---,Pp,/i))pi 



if d — p 



(3.19) 



\ Re{a p {p , pi, ■ ■ ■ , p 2 ,, p)) p P J 



if d = p+l. 



The following characterization of the mapping IT will be very useful later for computa- 
tional purposes: if G is an element of H^ +S (R K , T p ), then 



(UG)(p,p) = C- 1 -G( 1 ~ 1 -R,p), 



(3.20) 
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where 7 is any fixed element of T , R = (pi, p±, . . . , p p , p p ) if d = p and R = (p , Pi, Pi, ■ ■ ■ , P P , P P ) 
if d — p + 1, and where C is the following d x k matrix 



/ 1/2 1/2 
1/2 1/2 



C= { 



\ 







/ 1 

1/2 1/2 
1/2 1/2 



\ 



\ 


1/2 1/2 J 

• \ 

• 

• 

■ 1/2 1/2 J 



if d = p, k = 2p 



if d = p + l,K = 2p + l. 



(3.21) 



3.4 Parameter splitting 

There is a canonical direct sum decomposition of H^ +S {X) which will turn out to be quite 
useful for our purposes. Note that H^ +S {X) contains H^{X) (the //-independent polynomi- 
als) as a subspace, and consequently we can write 



H* +S (X) = H£(X) © P e K+s (X), 

where q G P £ K+S (X) if and only if q e H* +S (X) and q(x, 0) = 0. 

The homological operator C B (see (3.9)) preserves the decomposition (3.22): 



(3.22) 



C B (H*(R K )) C H*(R k ) and C B (P f K+s (R K )) C P f K+s {R K ). 



Moreover, 



H* +S (R K , T p ) = H£(R K , T p ) © P e K+s (R K , T p ), 



where H?(R K , T p ) = H*{R K ) n H* +S {R K , T p ) and P e K+s (R K , T p ) = P t 

If Cb\i and £b\2 represent respectively the restrictions of Cb on H^(R K ) and on P t 
then we have the following refinement of Proposition 3.5: 



(3.23) 

nm +s (R K ,T p ). 



K+sfrn>K\ 



Proposition 3.12 



H*(R K ) = H£(R K ,T p ) ©range C B \i, 
P e K+s (R K ) = P e K+s (R K , T p ) © range C B \ 2 - 
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The equivariant projection operator A defined in Definition 3.6 also preserves the decom- 
position (3.22), and we get the following refinement of Proposition 3.7 

Proposition 3.13 

A{H£{R K )) = H£(R K , T p ), A{P e K+s {R K )) = P e K+s {R K , T p ). (3.24) 
If A\i and A\ 2 represent respectively the restrictions of A on H^(R K ) and on P^ +S (R K ), then 

ker^lx = range Cb\i 

(3.25) 

kerA| 2 = range £5(2 



Remark 3.14 We note that there exist similar direct sum decompositions of H e +s (R d ) and 
of Hf +s (R d , Z 2)P ) using the subspace H d (R d ) C H°} +s (R d ) of ^.-independent polynomials 

H d+s (R d ) = H d (R d ) © P d+s (R d ) 

(3.26) 

Hf +s (R d , Z 2 , p ) = H d (R d , Z %p ) © P d + s (R d , Z 2 , p ) 

where q G P d+s {R d ) if and only if q G H d+s (R d ) and q(x,0) = 0, and where H d (R d ,Z 2 , p ) = 
Hf(R d ) n H d+s {R d ,Z 2 , p ) and P d+s (R d ,Z 2iP ) = P d+s (R d ) n H d+s (R d , Z 2 , p ) . Note that the 
mapping U defined in (3.18)-(3.20) preserves (3.23) and (3.26), i.e. 

Tl(H*(R K , T p )) = H d (R d , Z 2jP ) and U(P e K+s (R K , T p )) = P d+s (R d , Z 2jP ). 

Combining the results of this section with the Faria and Magalhaes normal form procedure 
described in section 2, we get the following version of Theorem 5.8 of [8] and Theorem 2.16 
of [9] which is adapted for our purposes 

Theorem 3.15 Consider the system (2.6) 

x = + fj(x,y,fi) 

d (3 ' 27) 
— y = ff(x,y,fj). 

i>2 

Write 

f}(x,Q, l i)=h j {x)+q j {x, l i), (3.28) 
21 



where hj e ifj"(M K ) and qj G P J K+S (M K ). TTien i/iere is a formal near-identity change of 
variables 

(x, y) — (£, y) + (C/ 1 ^), f/ 2 (£)) + (W\x, p), W 2 (x, p)) 

(where W l (x, 0) = 0, W 2 (x, 0) = which transforms (3.27) into system (2.8) (upon drop- 
ping the hats), and the flow on the invariant local center manifold y = is given by 

x = Bx + J2( (AUih, + YMx) + (A\ 2 (q 3 + Zj))(x : p)) (3.29) 
i>2 

where Y 2 = 0, Z 2 = 0, and for j > 3, Yj(x) and Zj(x,p) are the extra contributions to 
the terms of order j coming from the transformation of the lower order (< j) terms, and 
Zj(x, 0) = 0. 

4 Realizability: linear analysis 

In this section, we present the first of our main results on the realizability of the radial part 
(3.19) of toroidal normal forms (3.14) to any order for RFDEs (2.6) via the center-manifold 
normal form equations (3.29). 

We will define a linear operator between suitable spaces of polynomials, which arises in 
the context of the normal form transformations of (2.6). Our main result in this section will 
be to establish the surjectivity of this operator. Surjectivity will be the main ingredient in 
the proof of our main realizability results which will be presented in the next section. 

Again, we we will try as much as possible to use concise notation which will allow for the 
simultaneous treatment of both cases for A in Hypothesis 2.2(b). 

4.1 Preliminaries 

For given integers p > 1 and I > 2, and for k = 2p (respectively 2p+l) and d = p (respectively 
p+ 1), recall that Hf +s (R d , 1 2 , P ) is the linear space of homogeneous polynomials of degree 
t in the d + s variables p = (pi,...,p p ) (respectively p = (po, p±, . . . , p p ) = (po,p)) and 
p = (pi, . . . , p s ) with coefficients in M. d , and which are equivariant with respect to the Z 2jP 
action (3.17) on M d . Recall also that H^ +S (R K ) is the space of homogeneous polynomials of 
degree £ in the k+s variables ) (respectively x = (x , x\, x\, . . . , x p , x p )) 

and p = (pi, ...,p s ) with coefficients in R K , and H* +S (R K ,T P ) is the subset of H* +S (R K ) 
consisting of T p -equivariant mappings. 

Definition 4.1 Denote by V d+S (R) C H d+S (R) the subspace of homogeneous degree £ poly- 
nomials in the d + s variables v = (t>i, . . . , v p ) (respectively v = (v , v±, . . . , v p ) = (v , v) ) and 



22 



p = (f/,i, . . . , p s ) with real coefficients, spanned by the basis 
{ pflv 2k v c | c G {l,...,p}, (k, q) G N d+S , \q\ + 2\k\ + 1 = £} if al = p 

{ fi*v*> v 2k v c | c G {1, . . . ,p}, ((fco, k), q) G N d+S , \q\ + k + 2\k\ + 1 = 1} {J 
{ fi q VQ° v 2k I ((k , k), q) G Nq +s , \q\ + k + 2\k\ = £} if d = p + l, 

(4.1) 

where it is understood that if(k, q) = (ki, . . . , k p: qi, . . . , q s ) G Nq +s 7 £/ien /i 9 = {p\) qi • • • (p s ) 9s , 
■y2fc = -u 2fc = (t;!) 2 ^ 1 • • • (v p ) 2kp , \q\ = £V Qj an d \k\ = J2j fy- 

Note that, V d+S (R) is isomorphic to the vector space Hf +s (M, d : Z 2)P ), since this latter space 
has the following basis 

{ p q p 2k p c e c \ ce {l,...,p}, (k,q)eN d+s ,\q\ + 2\k\ + l=£} if d = p 

{ p q p k ° P 2k Pc e c+ i | c G {l,...,p}, ({ho, k),q)e N d+S , \q\ + k + 2\k\ + 1 = £} (J 
{p q p k °p 2k e 1 | ((A; ,A;),g)GNfMg| + A; + 2|A;|=£} if d = p+l, 

(4.2) 

where is a column vector with zeros on each row except the row, which is 1. Therefore, 
dim Hf +S (R d , Z 2)P ) = dim V d+S (R) . 

Since 5 is as in (3.2), then this corresponds to the following choice of basis for the center 
subspace P: 



*(t) 



(e iuJlt e- iuJlt ■ ■ ■ e*"** e - *^*) if d = p, k = 2p 

(1 e^ 1 ' e" iwi * ••• e^* e"^') if d = p + 1, k = 2p + 1. 
It follows that ^(0) in (2.5) is a k x 1 matrix 

{col(wi, «i, ...,-u p ,Mp") if k = 2p 
(4.3) 
col(u ,ui,ui, . . . ,u p , u p ) if n — 2p+l, 

where uq ^ is real and Uj ^ are complex, j = 1, . . . ,p. 
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4.2 Linear analysis 

Definition 4.2 Let S denote the normed real linear space of d x k matrices of the form 



( OL\\ OL\ 



l ' ' ' <Xi,p Qti,p 



\ 



«2,1 "2,1 • • • Ct2,p Oi 2 ,p 



if d = p, k = 2p 



M = < 



y o.p t i ctp t \ ■ ■ ■ ctp,p ®p,p J 



( 



a 0,0 "0,1 "0,1 ' ' ' a 0,p a 0,p ^ 

«i,o aTJ • • • 



(4.4) 



\ &pfl a p,l a p,l 



if d = p+ 1, K = 2p + 1, 



(where the a^o are rea/ and £/ie ajj, j > 1 are complex numbers) equipped with norm \ \M\ \ = 
max(]o;jj|). For any given M <E S , we define the ^-mapping associated to M 



J-t . TTd+S 
M ■ n t 



by 



(Ji I (h))(x,^ = ^(0)h(x(M) T ,fiy 
where x = (x±, x~i, . . . , x p , x^) (respectively x = (xo, xi,x±, . . . , x p ,x^) ). 
Let t — (ti, . . . ,Td) be a vector (as of yet unspecified) in M, d . Define 



E T = 



( \ 



r / 



v 



lLO\T p p — lLO\T v 



if d = p 



( 1 



e iu)pT! g-lWpTl \ 



{ V 



J 



if d — p + 1 , 
(4.5) 



and note that E T belongs to the space S of Definition 4.2. Define the linear mapping 
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by 

(£t(h))(x,n) = (Ji T )(h)(x,n) (4.6) 

where is the £- mapping associated to E T . 

Now, let n : H* +s (R k ,Jp) — ► Hf +s (R d ,Z 2 , P ) be the mapping defined in (3.18)-(3.20), 
and let A : H^ +S (R K ) — >■ H^ +S (R K , T p ) be the group averaging operator defined in Definition 
3.6. Our main result in this section is the following: 

Proposition 4.3 For an open and dense set U C R d , the following linear mapping is sur- 
jective for all r G U: 

IIoAo^: H d + S {R) — ► Hf +s (R d , Z 2jP ). 

Proof Let K be the d x d matrix such that Kj jk is equal to — 1 if j + k > d + 1 and is 
equal to 1 otherwise. It is easy to row reduce K to the identity matrix, so K is invertible. 
Therefore, K induces an automorphism of the space H d+S (R): 

K. : H d+S (R) — ► H d+S (R) 

(4.7) 

(JCh)(v^)=h(v (K) T ,fi). 

Define V e d+S (R) = JC' 1 (V e d+S (R)) , where V e d+S (R) is as in Definition 4.1, and let 

A/if : V e d+S (R) — > ^ +s (M d ,Z 2iP ) 

be the restriction of IT o Ao E^. to V/ +S (1R). Our approach to proving Proposition 4.3 will be 
to prove that there exists an open and dense set of points U C R d such that A/f is invertible 
for all t eU. 

If (A/if) is any matrix representation of A/if, then det((A/f )) is a real-analytic function of 
7i, . . . , Td (in fact, it is a polynomial in cos c^Tg and sin uj^q-, k £ {1, . . . g G {1, . . . , d}). 
Therefore, if we can show that det((A/if )) is not identically zero, the conclusion is a trivial 
consequence of this analyticity. This amounts to showing that there exists at least one point 
t* G R d such that with E T * as in (4.5), the mapping A/f, is invertible. We will prove this 
last claim with a sequence of five lemmas. 

Lemma 4.4 Let S be as in Definition 4-2. If M* G S is such that the restriction of the map 
UoAoJ^ to V d+S (R): 

UoAoJ^: V e d+S (R) — H d+s (R d , Z 2 , p ) 

is invertible, then there is a 5 = 8{M it ) > such that for all M in the 5-ball centered on M*, 
the restriction UoAoJ^: V e d+S (R) — ► Hf +s (R d , Z 2jP ) is invertible. 
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Proof This follows from the fact that the determinant of the map UoAoJ'j^ : Vf +S (R.) — > 
H d+s {R d , Z 2jP ) is continuous in the entries of M. I 

As mentioned above, \l/(0) in (4.3) is such that each of its components is non-zero. 
Therefore: 

Lemma 4.5 There exists a±, . . . , a p such that 

Re(e^ Uj )^0, j = l,...,p. 



Lemma 4.6 Let S be as defined in Definition 4-2 and a±, . . . , cr p be as in Lemma 4-5. Con- 
sider the following element I e S of the form (4-4) where ao,o = I in the case d = p + 1, 
and: 

( e i(y i if j = k > 1 

a i,k = { 

I if j^k 

If Jj : H d+S (R) — > Hp +S (R K ) is the £-mapping associated to I, then the restriction to 
V e d+S (R): 

UoAojf: V e d+S (R) — > Hf +s (R d , Z 2 , p ) 

is invertible. 



Proof We give only the proof in the case d — p, k — 2p. The other case (d — p + 1, 
k = 2p + 1) is treated in a completely similar manner. 

Consider the basis element p q v 2h v c of V[ +S (M) (see (4.1)). Then after an appropriate 
translation of the integration variables, we have 

(IloAojf)(ii*v 2k v c ) = 



(2tt)p Jo 



2tt 



2-k 



Q ( Pl (e- tei + e Wl )) 2kl ■ ■ • (p p (e-** + e^)) 2fe "(p c (e^ + e^))d9 1 ■ ■ ■ d9 p , 



-i0 c I a i0c\ 



where Q = C ■ diag(e i(?1 , e~ i(>1 , . . . , e* ', e~ ie ^) ■ diag(e i(Tl , e~ ieTl , . . . , e ia ", e~ ia ") • (0), (C as in 
(3.21) and \&(0) as in (4.3)). A simple computation then shows that 



(IloAojf)((i«v 2k - 



Re(e iac u c ) 



2k c + l 
k c +l 



(2h)\ (2k 2 )\ (2k p 



(k 2 \y (k p \y 
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fi q p 2k p c e c , 



where we remind the reader that e c is a p-dimensional column vector with zeros on each 
row except the c** 1 row, which is 1. Taking into account (4.1) and (4.2), under a suitable 
choice of bases for the spaces Vl +S (R) and H^ +S (R P , Z 2iP ), the matrix representation of the 
restriction of II o A o jf to V[ +S (R) is diagonal with non-zero diagonal entries. ■ 

Lemma 4.7 Let S be as in Definition 4-2, o~i, . . . , o~ p as in Lemma 4-5 and I as in Lemma 
4-6. Consider the element E* = KI G S, where the dx d matrix K is such that Kj^ is equal 
to -1 if j + k > d + 1 and is equal to 1 otherwise. If : H d+S (R) — ► H£ +S (R K ) is the 

£ -mapping associated to E m , then the restriction II o A o : Vf +S (R) — > H d+s (R d , Z 2jP ) is 
invertible. 

Proof Let / and jf be as in Lemma 4.6. Since E* = KI, it follows that jf^ = jf o JC, 
where K, is the automorphism defined in (4.7). So IloAoJ'^ = (IToAo j7/)o/C. Consequently, 

the restriction of II o A o to V e d+S (R) = IC' 1 (V e d+S (R)) is invertible. I 

Lemma 4.8 Let E* be as in Lemma 4-7 and let 5 = 5(E*) > be as in Lemma 4-4- 
There exists a r* G R p such that E T * in (4-5) satisfies \ \E T * — E*\\ < 5, and consequently if 
E** = 3ft t : Hf +S (R) — ► H^ +S (R K ) is the t-mapping associated to E T * , then the restriction 

Af^=UoAo £*, : Vf +S (R) — ► Hf +s (R d , Z 2;P ) is invertible (from Lemma 44). 
Proof Since the uj\, . . . , uj p are independent over the rationals, it follows that the set 

{(e^ 1 *, e iuJ2t , . . . , e*"**) \teR} 

is dense on the p-torus T p . Consequently, it is possible to choose a r* G R d such that each 
row of E T * is as close (in any given norm) as we wish to the corresponding row of E*. I 

The proof of Proposition 4.3 follows immediately from Lemmas 4.4-4.8. ■ 
4.3 Refinement 

In the next section, we will need a finer version of Proposition 4.3. If H d (R) denotes the 
subspace of /i-independent elements of Hf +S (R), then we have 

Hf +S (R) = H$(R) © P d+S (R), (4.8) 

where q G P £ d+S (R) if and only if q G H d+S (R) and q(v,0) = 0. Define Vf(R) = V e d+S (R) n 
Hf(R) and W e d (R) = V e d+S (R) n P d+S (R) and note that 

V e d+S (R) = V e d (R) © W? +S (R) (4.9) 
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is precisely the decomposition of Vf +s (Ei) into the direct sum of /i-independent elements of 
V e d+S (R) and elements of V e d+S (R) which vanish at fj, = 0. The automorphism K of Hf +s iR) 
defined in (4.7) preserves these decompositions, and we have 

V e d+S (R) = JC-\V^- a (R)) 

= /c- 1 (y/(R))e/c- 1 (vy/ +s (M)) (4.io) 

= V/(R) © W? +S (R) 

which is the decomposition of V e d+S (R) into the direct sum of ^-independent elements of 
V e d+S (R) and elements of V d+S (R) which vanish at fi = 0. 
Then, taking into account (3.26), we have 

Proposition 4.9 

dimt//(M) = dim^(M d ,Z 2iP ) and dim W d+S (R) = dimP d+s (R d , Z 2 , p ). (4.11) 
Furthermore, if r is as in Proposition 4-3, then 

(TloAo S e T )(V d (R)) = H d (R d , Z 2 , p ) and (n o A o S e T )(W d+s (R)) = P d+s (R d , Z 2 , p ). 



Proof We give the proof in the case d — p, the other case being treated in a similar 
manner. 

We note that 

\//(R) = span{t; 2 S c | c G {l,...,p}, k G N p ,2\k\ + 1 = 1}, 

and 

# P (M P ,Z 2 , P ) = span{p 2fe p c e c | cG {l,...,p}, k G Ng, 2\k\ + 1 = £}. 

Equation (4.11) follows from (4.9), (4.10) and the fact that V[(R) = /C _1 (V^ P (R)). It now 
follows from the theory presented in section 3 that 

(II o A o £ e T )(V e p (R)) C # P (M P , Z 2iP ) and (n o A o ^)(% +s (R)) C P P+S (M P , Z 2 , p ). 

The reverse inclusions then follow from the invertibility of jV^ and from (3.26), (4.8) and 
(4.11). ■ 
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5 Main Results 



We are now ready to state and prove our main readability results for both cases of Hypoth- 
esis 2.2(b), with the convention that respectively (d, k) = (p, 2p) and (d, k) — (p + 1, 2p + 1). 
It will be convenient to define the following linear spaces of (non-homogeneous) polynomials 

Definition 5.1 For an integer £ > 2, define 

V d+S (R) = ® l j=2 V d+S (R), V d (R) ee @] =2 V d (R), 

W d + S (R) ee ©J =2 W d+S (R) , H d £ +S (R d , Z 2jP ) ee ©J =2 H d+S (R d , Z 2iP ) , 

?#(R d , Z 2 , p ) ee ©j =2 tf/(R d , Z 2iP ), P £ d + s (R d , Z 2 , p ) ee ©j =2 P/ +s (R d , Z 2iP ), 

H d (R) ee ©j =2 if/(R), ee ©j =2 #f +s (R). 

Our first result addresses the issue of realizability of singularities and unfoldings within 
the class of scalar delay-differential equations with d delays. 

Theorem 5.2 Consider the RFDE (2.1), and let A denote the set of solutions of (2.3) with 
zero real part. Suppose that Hypothesis 2.2 is satisfied. Let t > 2 be a given integer. For 
each h E H d e (R d ,Z 2 , p ): 

h (p) = Yl Mp)> 

i=2 

(hj E H d (R d , Z 2)P ), j = 2, . . . , i) and each q E V d+s (R d , Z 2iP ): 

l 

Q(P,V>) = Y Qj(P,v), 
i=2 

(qj E P d+s (R d , Z 2)P ), j = 2, . . . , £), there are d distinct points Ti, . . . , r d E [— r, 0], an r) E 
V d (R): 

i 

viv) = Y ( 5 - x ) 

i=2 

fa E V d (R), j = 2, . . . ,£), and a £ E W d+S (R): 

£ 

av,fi) = YUv,P), (5.2) 
i=2 
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e W d+S (R), j = 2, . . . ,£), such that if 

F(z t , p) = r)(z(t + r 1 ),...,z(t + T d )) + £(z(t + r 1 ),...,z(t + T d ),p) 

in (2.2), then in polar coordinates, the radial part of the center manifold equations (3.29) in 
T p -equivariant normal form up to degree t reduces to p = h(p)+q(p, p), where p = (pi, . . . , p p ) 
or p = (p , pi, . . . , p p ). In fact, r can be chosen in an open and dense set of [— r, 0} d , 
independently of the particular h and q to be realized (i.e. only i] and £ must be changed in 
order to account for different jets to be realized). 

Proof Choose a point r G [— r, 0} d such that the previously defined linear mappings 

Mi : V d+S (R) — > H d+s (R d ,Z 2 , p ) 

are invertible for all j = 2,...,£ (from Proposition 4.3, this is possible for an open and 
dense set of points in [— r, 0} d ). Suppose r\ is an arbitrary polynomial of the form (5.1), £ is 
and arbitrary polynomial of the form (5.2), and suppose F in (2.2) is such that F(z t ,p) = 
rj(z(t + Ti), . . . , z(t + Td)) + £(z(t + Ti), . . . , z(t + Td), p). Using Theorem 3.15, it is possible 
to define successively at each order near identity changes of variables of the form 

(x,y) = (x,y) + (Uj(x) +W}(x,p),Uf(x) +Wf{x,p)), (5.3) 

where Wj(x,0) = 0, % = 1,2, which transform (2.6) into (2.8), and the center manifold 
equations are as in (2.9), with 

g\(x,Q,v) = A(£* V2 )(x) + A(£^ 2 )(x,p) 
gl(x,0,p) = A (£3 rj 3 + Y 3 ) (x) + A £ 3 + Z 3 ) (x, p) 

; (5.4) 

g](x,0,p) = A(S^ 3 +Y 3 )(x) + A(SU, + Z 3 )(x,p) 

In (5.4), Si are as in (4.6), A is the T p averaging operator (3.6), and Yj(x) and Zj(x,p) are 
the extra contributions to the terms of order j coming from the lower order (< j) changes 
of variables, and Zj(x, 0) = 0. Hence, the terms Yj and Zj are completely determined once 
the normalizing procedure arrives at order j. More precisely, Yj is determined explicitly 
in terms of r] 2 , . . . , rjj-i, U 2 , . . . , Uj_ ± , i — 1,2 and Zj is determined explicitly in terms of 
T)2, ... , r)j-i, £2, • • • , U 2 , . . . , Uj_ ± , W 2 , . . . , Wj_i, i — 1,2. Taking into account (5.4) and 
using the convention Y 2 = 0, Z 2 = 0, the center manifold equations (2.9) are TP-equivariant, 
and in polar coordinates, the uncoupled radial part (truncated at order £) is of the form 

l 

p = E W n + ( n a )( y Mp) + W Si + (n o a)(zMp, p)} ■ 

i=2 
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Thus, using Proposition 4.9, we get the desired result if we set 

Vj = (A/?)' 1 (hi - (n o AM)), = (A^)' 1 (qj - (H o A)^)), j = 2, . . . , £. (5.5) 

■ 

Theorem 5.2 has an important interpretation in terms of the singularity and unfolding 
theory of scalar delay-differential equations. Suppose (2.1) satisfies the hypotheses of The- 
orem 5.2. Let h(p) be any given (parameter independent) element of 7if (M. d , Z 2iP ), ^ > 2. 
Then Theorem 5.2 implies that (under generic conditions on T\, . . . , T d ) there exists an un- 
parameterized nonlinear polynomial delay-differential equation 

z(t) = L z t + r]{z{t + n),...,z(t + r d )) (5.6) 

whose dynamics on a center manifold up to order £ have as uncoupled radial equations 
p = h(p). Therefore, generically, any finitely- determined singularity within the space of 
Z 2;P - equivariant radial equations can be realized by an appropriate choice of rj in (5.6). 

Now, suppose that h G H d+s (M. d ,Z 2 , p ) is an equivariant unfolding of the finitely-determined 
singularity h above, i.e. h(p, p) is such that h(p, 0) = h(p). Then q(p, p) = h(p, p) — h(p) is an 
element of Vf +s (M, d , Z 2)P ). Theorem 5.2 implies that there exists a parameterized nonlinear 
polynomial delay-differential equation of the form 

z t = L z t + rj(z(t + Ti), . . . , z(t + r d )) + £(z(t + n),...,z(t + T d ),n), (5.7) 

with £(z(t + Ti), . . . , z(t + T d ), 0) = 0, whose dynamics on a center manifold up to order I 
have as uncoupled radial equations p = h(p,p) = h(p) + q(p,p). Therefore, the unfolding 
r)(z(t + Ti), . . . , z(t + Td)) + £(z(t + ri), . . . , z(t + r d ), p) of i] realizes the unfolding h(-, p) of 
the singularity h on the center manifold. 

In the theory of classification of singularities of equivariant vector fields [12, 13], one 
often defines a suitable equivalence relation on a given space of vector fields (by requiring 
preservation of certain local qualitative features of the flow associated to the vector field), 
and then classifies the equivalence classes in terms of a (hopefully finite) set of conditions of 
the Taylor coefficients of the vector field. One then wishes to characterize the "likelihood" 
of a given singularity, /, by computing its codimension, which roughly speaking, is the codi- 
mension of the equivalence orbit through /. Finally, one then uses this idea of codimension 
to construct a versal unfolding of the singularity (perturbing in transversal directions to the 
equivalence orbit). 

Suppose / : R d — > R d is a smooth vector field vanishing at the origin and equivariant 
with respect to the group Z 2 p previously defined. Typically, the computation of codimension 
of the singularity / is done by first identifying a polynomial h G Hf (R d , Z 2>p ) (for some 



31 



suitable £) which is equivalent to / regardless of the Taylor coefficients of / of order greater 
than I. Then, one constructs the tangent space within 7i d (R d , Z 2j? ) to the equivalence 
orbit of h through h, T h C ^(R rf , Z 2jP ). Finally one finds a complementary subspace 
C h C H d e (R d , Z 2jP ) such that H d e (R d , Z 2 , p ) = T h © C h . The codimension of / is then identified 
with the dimension of Ch (i.e. codim/ = codimT/j = dimC^). We say that / is generic with 
respect to the equivalence relation if the codimension of / is zero. The following theorem 
addresses this issue within the context of realizability. 

Theorem 5.3 Consider the RFDE (2.1) in the unparametrized (s = 0) case, and let A 
denote the set of solutions of (2.3) with zero real part. Suppose that Hypothesis 2.2 is satisfied. 
Suppose that the nonlinear term F(z t ) is of the general form 

F(z t ) =r ] (z(t + T 1 ),...,z(t + T d )), 

where 77 is smooth. Then the local dynamics of (2.1) near the origin on an invariant center 
manifold can be described by a system of ordinary differential equations on R K . Moreover, 
this ODE system can be brought into T v -equivariant normal form to any desired order t, and 
the resulting (truncated at order t) normal form can be uncoupled into two sub-systems 

P = h(p;V, r ) (5-8) 

= k(p; V ,r), (5.9) 

where r = (n, . . . , r d ) e R d , h{- ; 77, r) G H d e {R d , 1 2 , P ) and k{- ;r],T):R p — ► R p . For given 
r G R d , consider the following mapping: 

F T :H d e (R) — ► H d (R d ,Z 2 , P ) 

r] 1 — ► Ft(v) = K;v,t), 

where h is as in (5.8). Then there is an open and dense setU C R d , such that for all r eW, 
T T is a submersion. Consequently, if M. C 7Y^(lR d , Z 2)P ) is a smooth manifold, then for all 
r GW, JF j t 1 (7W) is a smooth submanifold ofHf(R), and 

codim^ 1 ^) = codimA^ (5.10) 



Proof The fact that the center manifold equations are given by (5.8) and (5.9) has already 
been proved. 
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The mapping T T is computable similarly to (5.5): if r\ — Y^ 3 =2 Vji w ^ n Vj e H d (R), then 

i 

tm = K- ■ v , t) = ((n o a o £i)( Vj ) + (n a){y 3 )\ 

J'=2 

where Y 2 = and is a smooth function of r) 2 , . ■ ■ ,r]j-i f° r 3 > 2. Thus, if C — X^=2 G> 
with 6 G Hf(M), then 



^r(r)) • C = E ( ( n ° A ° £ r)(0) + o A) (j2 Y Md 

j=2 \ \i=2 

where Yji — if % > j. From Proposition 4.3, there is is an open and dense set W C l d 
such that for all r G U, DT T {r]) is onto Tif (M. d , Z 2jP ), and consequently T T is a submersion. 
Equation (5.10) follows from the transversal mapping theorem [1]. I 

The next result states that the number of delays, d, shown above to be sufficient to realize 
any arbitrary element of TifiW 1 , %2,p), is optimal for that purpose. 

Theorem 5.4 Consider the RFDE (2.1) in the unparameterized (s = 0) case, and let A 
denote the set of solutions of (2.3) with zero real part. Suppose that Hypothesis 2.2 is satisfied. 
Suppose that the nonlinear term F(z t ) is of the general form 

F(zt) = ri(z(t + r 1 ),...,z(t + r d _i)), 

where rj is smooth. Then the local dynamics of (2.1) near the origin on an invariant center 
manifold can be described by a system of ordinary differential equations on M K . Moreover, 
this ODE system can be brought into T p -equivariant normal form to any desired order I, and 
the resulting (truncated at order t) normal form can be uncoupled into two sub-systems 

p = Hp;v, t ) (5-11) 

6 = k(p; V ,r), (5.12) 

where r = (n, . . . , r d _i) G R d -\ h(- ; r], r) G H d (R d , Z 2 , P ) and k{- ; rj, r) : W — ► W. For 
given r G consider the following mapping: 



Fr-.H^iR) — > H d e (R d ,Z 



T) i — > J- r (rj) — h(y- ; tj,t), 



(5.13) 



where h is as in (5.11). Then there is an integer £ > 2 such that T T is not surjective if 
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Proof It will be sufficient to show that for fixed d, 

dimHt 1 



dimH d (R d , Z 2 , 



0(r v ) (5.14) 



as I — > oo. First, note that it is well-known [14] that for given integers m > 1 and £ > 1, 
the number of solutions in non-negative integers for the equation 

&i -|- • • • -\- k m — I 



is 



m + £-l 
m — 1 



Thus, 



dim^-^R) = dimHf-\R) = ^( rf ^ 2 2 ) 

i=2 i=2 ^ ' , 

(5.15) 

= (^ d ~ d ]_Y^) -d = 0{£ d ' 1 ) as £^00. 
Using a similar (but slightly lengthier) computation, we can show that 

dimH d (R d , Z 2)P ) = 0(f*) as £ -> 00, 

which establishes (5.14) and concludes the proof of this theorem. I 

Theorem 5.4 is important in the problem of establishing whether or not there are restric- 
tions on the possible phase portraits for an unfolding of a given singularity h G Hg(M. d , Z 2;P ) 
when such an unfolding arises from center manifold reduction and phase/amplitude decou- 
pling of a nonlinear delay-differential equation (2.1). This theorem allows one to conclude 
that, at least for £ large enough, such restrictions are likely to occur if the number of delays 
in the nonlinearity F in (2.2) is less than d. For example, this question was addressed in [4] 
in the context of the non-resonant double Hopf bifurcation. In this case, they show that if 
the nonlinear part of (2.1) contains 2 delays, then generically any cubic order radial equation 
(5.8) can be realized by appropriate choice of the nonlinear coefficients in (2.1) (note that 
our Theorem 5.2 recovers and generalizes that result). However, in [4], it is also shown that 
if the nonlinear part of (2.1) depends on only one delay, then not all equivalence classes 
of phase portraits in the versal unfolding of the radial equations (5.11) can be attained by 
variation of the nonlinear coefficients in (2.1), for fixed values of oui, uo 2 and r. The next 
example treats this specific case in the context of Theorem 5.4. 
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Example 5.5 In the case where A = {±iui, ±icu 2 }, we get from (5.15) that 

dim^(R) =£-1. 

It is also easy to show that if £ = 2L + j, where L > is an integer and j G {0, 1}, then 

dimft 2 (R 2 ,Z 2 , 2 ) = L(L + 3). 

Thus, from I — 3 onward, we have dim Ti\ (R) < dim7i 2 (R 2 , Z 2)2 ). In particular, dim 7^3 (R) = 
2 and dimW§(R 2 , Z 2)2 ) = 4. Therefore, at cubic order, the mapping jF r in (5.13) 

F T :Hl(R) — > ^i(R 2 ,Z 2 , 2 ) 

77 1 — ► ^(77) = ; 77, r), 

is not surjective, and so there are elements of 7i§(R 2 , Z 2j2 ) which can not be realized by 
any element of Ji\ (R). In fact, T T (7£j(R)) is a two-dimensional smooth surface in the 4- 
dimensional space 7i|(R 2 ,Z 2i2 ), such that •Tv(O) = 0. Specifically, if we write the general 
element of Hl(M) as 

b 2 v 2 + b 3 v 3 
and the general element of 7i 2 (R 2 ,Z 2j2 ) as 

(a u pj + aupQpi 
{a 2 ipj + a 2 2p 2 2 )p2, 

then the mapping T T can be represented by the following mapping from R 2 into R 4 : 

aij(b 2 ,b 3 ) = aijb 2 2 + pijb 3 , ij = 1,2 (5.16) 

where the real coefficients and are determined from r, uo\ and uo 2 . 

Note however that the problem of determining whether or not there are restrictions is 
somewhat more subtle than one of surjectivity, since the topological types of the possible 
phase diagrams in the unfolding space for the double Hopf bifurcation are determined by 
the sign of the cubic coefficients in the radial equations (and not their actual values). In 
the (b 2 ,b 3 ) plane, the zero level sets of the in (5.16) are (at most) four distinct curves 
(parabolae generically) which intersect only at the origin. Consequently, there are at most 
four distinct open regions in the (b 2 ,b 3 ) plane in which the signs of the coefficients are 
constant and non-zero. It is then easy to see that it is impossible to realize the twelve 
possible sign combinations (see [15]) which characterize the complete unfolding space of the 
double Hopf bifurcation, and so there will be restrictions on the phase portraits when the 
nonlinear terms in (2.1) contain only one delay. 
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For general RFDEs (i.e. not necessarily delay-differential equations), we have the follow- 
ing result on realization of unfoldings: 

Theorem 5.6 Consider the general nonlinear RFDE 

z{t) = L z t + N(z t ) (5.17) 

where L : C — > R is a bounded linear operator from C = C ([— r, 0] , R) into R, and iV zs 
a smooth function from C into R, witt iV(0) = 0, DN(0) = 0. Le£ Ao denote the set of 
solutions of (2.3) with zero real part and suppose that Hypothesis 2.2 is satisfied. Then the 
local dynamics of (5.17) near the origin on an invariant center manifold can be described 
by a system of ordinary differential equations on R K . Moreover, this ODE system can be 
brought into T p -equivariant normal form to any desired order I, and the resulting (truncated 
at order £) normal form can be uncoupled into an uncoupled d-dimensional system and a 
p-dimensional system 

p = h(p;N) (5.18) 
9 = k(p;N), (5.19) 

where for given N, h(- ; N) is some element of Hf(M. d ,Z 2 , P ), and k(- ; N) : W — >W. 
Let h(p, p) be an s-parameter equivariant unfolding of h of degree at most i, i.e. h e 
H.e +s (M. d , Z 2jP ) and h(-,0) = h(- ; N). Then there exists an s-parameter unfolding of (5.17) 
of the form 

z(t) = L (zt) + N(z t ) + £(z(t + n),...,z(t + T d ),p) (5.20) 

(where r = (ri, . . . , r<j) G R d ; and I; e W^ +S (R) vanishes at p = 0) which realizes the unfolded 
radial equations 

p = ~h{p,p) 
on an invariant center manifold for (5.20). 

Proof Choosing Wj = 0, Wj = and choosing Uj and Uj appropriately in (5.3), the 
center manifold equations for (5.17) truncated at order £ are equivalent to (5.18) and (5.19). 

Using Theorem 3.15, for arbitrary £ = X^=2 ^ VV/ +S (R), there is a sequence of near- 
identity changes of variables (5.3) (with Uj and Uj as above) for which the uncoupled radial 
part of the center manifold equations for (5.20) truncated at order £ are 

i 

p = h(p;N) + J2M^ + (noA)(Z J )) 

j=2 
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where (II o A)(Zj) is some known element of P d+s (R d , Z 2iP ). The conclusion follows from 
setting 

^ = {^Y\h-h-{YioA){Z 3 )). 



6 The ±iu;, (dz iu>i, dz 20*2) and (0, ± iu;) Singularities 

Our results in Theorems 5.2, 5.3 and 5.6 allow us to recover some previous results on realiz- 
ability and (lack of) restrictions for Hopf bifurcation, non-resonant double Hopf bifurcation, 
and the (0, ±iu) singularity in scalar RFDEs [11, 4]. 

Corollary 6.1 (Theorem 1 of [11]) Consider the RFDE (2.1) in the unparameterized case 

z(t) = L z t + F(z t ), (6.1) 

such that the characteristic equation (2.3) has simple purely imaginary roots ±iuj ^ and 
no other roots on the imaginary axis (simple Hopf bifurcation) . If 

F(z t ) = A 2 (z(t + r)f + A 3 (z(t + r)) 3 , r G [-r, 0] (6.2) 

then the uncoupled radial part of the center manifold equations to cubic order are 

p = ap 3 , (6.3) 

where a = a(A 2 , A 3 ; r, uj). Generically, the non- degeneracy condition a ^ is satisfied. In 
fact, for any a e 1, (6.3) can be realized with A 2 = for an appropriate choice of A 3 in 
(6.2). Furthermore, in the case a ^ 0, the versal unfolding 

p = pp + ap 3 

of ( 6. 3) is generically realized ( modulo a rescaling of the parameter) by the following unfolding 
of (6.1) 

z(t) = L z t + F(z t ) +pz(t + r). 

Corollary 6.2 (Theorem 3.1(1) of [4]) Consider the RFDE (2.1) in the unparameter- 
ized case 

z(t) = L z t + F(z t ), (6.4) 
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such that the characteristic equation (2.3) has simple non-resonant purely imaginary roots 
±iu>i, ±iu>2, and no other roots on the imaginary axis (non-resonant double Hopf bifurcation). 
U 

F(z t ) = A 20 (z(t + T l )) 2 + A 11 z(t + T 1 )z(t + T 2 ) + A 02 (z(t + T 2 )) 2 + 

A 30 (z(t + n)) 3 + A 21 (z(t + n)) 2 z(t + r 2 ) + A 12 z(t + n)(z(t + r 2 )) 2 + (6.5) 
A 03 (z(t + T 2 )f, 

where t±, t 2 G [— r, 0], then the uncoupled radial part of the center manifold equations to cubic 
order are 

Pi = (a n p 2 i + a 12 pl)pi 

(6.6) 

p 2 = {a 2l p\ + a 22 pl)p 2 , 

where = aij(A 20 , A n , A 02 , A 30 , A 2 i, A 12 , A 03 ; t±, t 2 , u>i,u 2 ). Generically, the non- degeneracy 
condition a\\a 22 — a 2 \a\ 2 ^ is satisfied. In fact, for any an, di 2 , a 2 i, a 22 e M, (6.6) can be 
realized with A 20 = An = A 02 = for an appropriate choice of A 30 , A 2 i, Ai 2 , A 03 in (6.5). 
Furthermore, in the case a\\a 22 — a\ 2 a 2 \ ^ 0, the versal unfolding 

Pi = (pi + a u pj + a 12 p 2 2 )p 1 

(6.7) 

p2 = (P2 + a 21 pl + a 22 pl)p 2 

of (6.6) is generically realized (modulo a linear change of parameters) by the following un- 
folding of (64) 

z(t) = L z t + F(z t ) + nx z(t + Ti) + p 2 z(t + T 2 ). 



Remark 6.3 We would like to clarify the statement "modulo a linear change of parameters" 
in the preceding Corollary. According to the notation we have established in this paper, we 
have 

Vi +2 (R) = V 2 2+2 (R) = W 2+2 (R) = W 2+2 (R) = spanf//!^, p x v 2 , p 2 v u p 2 v 2 } 

and 

H 2 2 + 2 (R 2 , Z 2 , 2 ) = P|+ 2 (R 2 , Z 2 , 2 ) = span [to ^ ) , to ( ° ) ,to ( * ) , to ( ° ) } ■ 
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From Proposition 4-3, the mapping 



TV 2 : V 2 2+2 {R) — ► H 2+2 {R 2 ,Z 



2,2 J 



is generically invertible. Since the mapping Af 2 does not have any effect on the parameters 
/ii and /i 2; generically we have 

(AC 2 )' 1 [pj ( ^ )) = to ( muVl + m i2^2), j = 1,2 

(Af 2 )- 1 [pj ( p 2 = Pi (mzivi + m 22 v 2 ), j = 1,2, 

and it follows that the 2x2 matrix M = (my) invertible. Consequently, the unfolding 
(6.7) of (6.6) is realized by the following unfolding of (6.4) 

z(t) = L z t + F(z t ) + pi (m n z(t + n) + m 12 z(t + r 2 )) + /i 2 (m 2 iz(t + n) + m 22 2:(t + r 2 )). 
W^e oei £/ie conclusion of Corollary 6.2 by performing the linear change of parameters 

fix = mufXi + m 21 fi 2 , Ji 2 = m 12 pi + m 22 fi 2 

and dropping the tildes. 

Corollary 6.4 (Theorem 2 of [11]) Consider the RFDE (2.1) in the unparameterized case 

z(t) = L z t + F{z t ), (6.8) 

such that the characteristic equation (2.3) has simple purely imaginary roots ±iu ^ 0, a 
simple root at 0, and no other roots on the imaginary axis (interaction of a simple bifurcation 
and a Hopf bifurcation). If 

F(z t ) = A 20 (z(t + n)) 2 + A n z(t + n)z(t + r 2 ) + A 02 (z(t + r 2 )) 2 (6.9) 

where r 1? r 2 e [— r, 0], then the uncoupled radial part of the center manifold equations to 
quadratic order are 

Po = hpl + hp 2 

(6.10) 

p\ = ctipoPi, 

where the coefficients a±, b\ andb 2 are functions of (A 20 , An, A 02 ; t±, t 2 , u). Generically, the 
non- degeneracy conditions a\ 7^ 0, b\ 7^ 0, b 2 7^ and a\ 7^ b 2 are satisfied. 
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Corollary 6.5 Consider the singularity (6.10) in the non- degenerate case a\ ^ 0, bi ^ 0, 
b 2 7^ and a x ^ b 2 . Then the following Langford unfolding [21] of (6.10) in the transcritical 
case 

Po = PiPo + hpl + b 2 p\ 

(6.11) 

Pi = P2P1 + OiPoPi, 

is generically realized (modulo a linear change of parameters) by the following unfolding of 
(6.8) 

z(t) = L z t + F(z t ) + /Xi z(t + Ti) + fl 2 z(t + T 2 ). 

7 Conclusions 

We have established a framework for the realizability problem for scalar RFDEs which ex- 
ploits fully the toroidal equivariance of normal forms of bifurcations associated with purely 
imaginary eigenvalues. This has allowed us to recover and significantly generalize recent 
results of Faria and Magalhaes [11] and of Buono and Belair [4]. As mentioned in the Intro- 
duction, it is important for modelers using RFDEs to be able to accurately assess the range 
of possible dynamics accessible within their models. For this purpose, this paper gives a thor- 
ough analysis of this question in the case where the model is a nonlinear delay-differential 
equation undergoing non- resonant multiple Hopf bifurcation or transcritical/non-resonant 
multiple Hopf interaction. Specifically, we split the dynamics of the normal form into com- 
ponents which are normal to the orbits of a torus group, and components which are tangent 
to these group orbits. Sharp estimates on the number of delays are then given for the re- 
alizability of the normal "radial" part of the normal form by nonlinear delay-differential 
equations. The case of saddle- node/non- resonant multiple Hopf interaction will be treated 
using similar techniques in a subsequent paper [5]. 

The generalizations we have achieved in our paper are twofold. First, we can treat within 
a unified framework the general case of p non-resonant Hopf eigenvalues and the interaction 
between simple steady-state bifurcation and p non-resonant Hopf bifurcation. Second, in 
contrast to [11] and [4] where only the generic (non-degenerate) cases are treated, we can 
treat the general finitely-determined case (whether degenerate or not) and its unfoldings, 
also within a unified framework. Note that in parameterized families of vector fields with 
sufficiently many parameters, it becomes possible to violate any specified non- degeneracy 
condition which is expressed in terms of the Taylor coefficients of the vector field up to 
some finite order. Therefore, it becomes desirable to have a framework in which these 
degenerate cases and their unfoldings can be systematically treated. Our results provide 
such a framework. 
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Open problems of interest related to this analysis and worthy of further investigation are 

• relaxing the restriction to scalar RFDEs in order to consider n > 1 dimensional systems 
of RFDEs 

• incorporating resonances in the purely imaginary eigenvalues and repeated eigenvalues 
with Jordan blocks. 
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A Proof of Proposition 3.5 

Let / be a given element of H* +S (R K ), and consider / = (/, 0) G H^ +S (R K+S ). From Propo- 
sition 3.2, there exists h = (h,h 2 ) G H* +S (R K+S ) and a unique g = (gi,g 2 ) G H* +S (R K+S ,T) 
such that 

(f,0) = C 6 (h 1 ,h 2 ) + (91,92). (A.l) 

Now, C B (hi,h 2 ) = (Cshi, D x h 2 Bx), so it follows that g 2 = —D x h 2 Bx. Consequently, (A.l) 
can be rewritten as 

(/,0) = (£ B /n,0) + (<7i,0), 

and thus 

f = C B h 1 +g 1 , (A.2) 

where g x G H* +S (R K , T p ). So H« +S (R K ) = H« +S (R K , T p ) +range£ B . Suppose / = in (A.2), 
then it is easy to see that £§(/ii,0) + (<7i,0) = (0,0), and from Proposition 3.2, it follows 
that gi = and C B h\ = 0. Therefore, 

H* +S (R K ) = H« +S (R K , T p ) © range C B 
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B Proof of Proposition 3.7 

For a given / G H* +S (R K ), let g = Af; then 

(Ag)(x,n) = / 7^(7" 1 x,/i)rf7= / 7 (/ 7 /(7 _1 7 -1 a;, //) cfy j 
</r -/r \^r / 

= / 77 /((tt)" 1 ^/^)^ ^7 

= / ( / lf(l~ 1 x,^)dl ) d^i = / 7 f(^ 1 x,fi)d- f 
Jr \Jr J Jr () 

= (Af)(x,n), 

where the second to last line holds because of the translation invariance and the normalization 
of the Haar integral. So A is a projection. 
Now, let / G ranged, then Af = /, i.e. 

f(x,fj) = / 7/(7~ 1 a;,^)rf7. 
</r 

So, for any a G r , we have 

af{a~ 1 x^) = a / 7/(7"V -1 x, //) cfy = / ^7 /((<77) _1 a;, /i) cfy 
•/r -/r 

= / 7 f(l' 1 x,^)d-f = f(x). 

Therefore, / G H* +S (R K , T p ). On the other hand, if / G #; +s (M K , T p ), then 

(Af)(x,fi)= 7 /( 7 - 1 x,/i)rf7= / f(x,fi)d'y = f(x,(i), 
J t Jr 

so / G range A. This establishes (3.10). We now establish (3.11). Since A is a projection, 
then 

H* +S (W) = range A © ker A. 

From Proposition 3.5, we conclude that dim ker A = dim range Cb- Thus, we need only 
show that range Cb C ker A. In order to show this, we will need the following 

Lemma B.l Let g : T — ► M. K be a continuous function, then 



f g{ 1 )d 1 = lim 1 I g(e Bs )ds. 
Jr T ^°° 1 Jo 
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h(9i, . . . ,9 P ) 



Proof For a given q G {1,2,... ,p}, consider the rotation matrix 

Rg{0) = diag{e l ^ e ,e~ iuJ " e ) 

which is T q = 27r/cJ 9 -periodic in 9. Then R q (9) = R q (9T q / '(2tt)) is 27r-periodic in 9. By 
hypothesis, the set 

J 2tt 2tt1 

til'"'?;/ 

is algebraically independent. Let T' p = [0, 27r] p , then we can parameterize T as follows: 

/i : T p — > T 
diag^i^i),...,^^)) if « = 2p 

diag(l, J R 1 (^ 1 ),..., J R p (^)) if « = 2p+l. 

Define g : T p — > R K by . . . , 9 P ) = (go h)(9 1: . . . , 9 P ). Obviously, g is 27r-periodic in 
each of its entries, and 

— — / •••/ ~g(9 1 ,...,9 p )d9 1 ---d9 p = g( 1 )d 1 . 
Noting that h . . . , ^) = e Bs and using Lemma 4.1, P. 430 of [6], we get that 

" $SL¥ Jo g(eB ' )ds ' 

which yields the desired result. ■ 
Now, let / G range then there exists g G i7" +s (R K ) such that 

D x g(x, n)Bx - Bg(x, /i) = f(x, (i), V (x, n) G K K+S . 
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Therefore, using Lemma B.l, we get 



{Af)(x,n)= [ 7/(7^^)^7= lim i [ T e Bs f(e- Bs x,n)ds 
Jr T ^°° 1 Jo 

1 /" T 

= lim - / e Bs (D x g(e- Bs x,ii)Be- Bs x - Bg(e~ Bs x, fi)) ds 

T^oo T J 

= l im I f T A. (e Bs g(e- Bs x ,//)) rfs 

e BT g(e~ BT x,n) -g(x,n) 
/"•x " T 

and this last limit is equal to 0, since the numerator is bounded in T for any given (x, /i) G 
R K+S . So we conclude that / G ker A, and thus that ker A = range C B - This establishes 
(3.11), and concludes the proof of Proposition 3.7. I 
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